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Tl' ABSTRACT 

^ : 

■ Geometrically thick disks may form after tidal disruption events, and rapid 
^ ■ accretion may lead to short flares followed by long-term, lower-level emission. 

Using a novel accretion disk code which relies primarily on global conservation 
^ ■ laws and the assumption that viscosity is everywhere positive, a broad range 

■ of physically allowed evolutionary sequences of thick disks is investigated. The 
CN ' main result is that accretion in the thick disk phase can consume only a fraction 
Q ■ of the initial disk material before the disk cools and becomes thin. This fraction 

■ is ~ 0.5 — 0.9 for disruptions around 10^ to lO^M© black holes and is sensitive 
. to the mean angular momentum of the disk. The residual material will accrete 

Q_|. in some form of thin disk over a longer period of time. The initial thick disk 

Q . phase may reduce the dimming timescale of the disk by a factor of ~ 2 from 

^ . estimates based on thin disks alone. Assuming an 0.5 Mq initial thick disk, 

. even if the thin disks become advection dominated, the black hole mass to light 

>• . ratio can rise above Mq/ Lq = 1 in no less than 20 (0.1 /a) to 2000 (0.1 /a) years 

^ . following a tidal disruption event, depending on the mass of the black hole and 

the initial conditions of the encounter. The long-term emission will be most 
prevalent around lower mass, IO^Mq black holes. If the tidal disruption rates 
in these galactic nuclei are ~ 10~^ yr~^, then about 10% of the nuclei should 
exhibit the long-term UV/optical emission at a level of ~ 10^^ ergs s~^. 
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1. Introduction 

If black holes reside in the centers of galaxies, one clinching sign of their existence would 
be the tidal disruption of stars when a star passes within a tidal radius, Rt ~ 25Mg~^^^i?g, 
of the black hole, where Mq is the black hole mass in units of IO^Mq and -Rs is the 
Schwarzschild radius. Tidal disruption events are rare (~ 10~^ per year per galaxy), 
bright (L ~ -^^Edd), short (months to years), and seemingly unavoidable consequences of 
10^ — IO^Mq black holes in galactic centers. The general picture for tidal disruption of a 
star has been developed over the last 20 years (e.g. Hills 1975; Young, Shields, & Wheeler 
1977; Rees 1988). The picture has evolved from one in which the disruption rate was so high 
as to fuel quasars to one with a more modest rate in which disruptions are rare transients 
which may nevertheless, slowly grow a lower mass, ~ 10'' M0, black hole (Goodman & Lee 
1989). 

At the time of disruption, half of the star is captured onto a range of elliptical orbits 
around the black hole; the remainder is ejected from the system. Due to the relativistic 
effect of orbital precession, debris orbits cross and may form strong shocks ( |Rees 1988]) . 



The explicit evolution of the debris is complex and is dependent on the spin and mass 
of the black hole as well as the properties of the star and the stellar orbit. As discussed 
by Kochanek (1994), at intersections between the debris streams, energy and angular 
momentum may be transfered. The computational complexities of the problem prohibit a 
direct numerical solution, but as an approximation, a large fraction of the debris can be 
said to circularize in a time, 

^cir ^orb^min) (1) 

where number greater than 1, and probably between 2 and 10. The minimum 

return time, tmin is the time for the most bound material to return to pericenter, Rp is 



where i?^ and are the radius and mass of the disrupted star. 

The circularization process may produce a thick torus which would accrete on a 
fast timescale. Initially, the disk is expected to be geometrically thick and radiation 
pressure supported, because of the large energy difference (~ GMM^/8Rp) between the 
nearly unbound elliptical orbits and a subsequent circular orbit at ~ 2i?p, which is set by 
conservation of angular momentum. The disk is also highly optically thick with r > 10^. In 
comparison, our sun has r ~ 10^^ and a radius ~ 100 times smaller than the characteristic 
radius of a thick disk. The thick disk's structure is very much like the envelopes of 
massive stars, except that self gravity plays no role; the gravity in the disks is dominated 
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by the black hole. The disk will remain geometrically thick if energy is dissipated fast 
enough. In general, this energy dissipation rate corresponds to an accretion rate which is 
super-Eddington. The physical parameters of an tidal disruption encounter which produce 
super-Eddington accretion are described in Ulmer (1997; hereafter U97). Delineating the 
region in which thick accretion occurs is important not only because super-Eddington tidal 
disruption events are the brightest and therefore, the most readily observed flares, but also 
because one expects the accretion of the stellar debris to be fastest when the disk is very 
thick. In thin disk theory, the accretion time scales as a~'^{z/r)~'^ , where a is the viscosity 
parameter, z is the height of the disk, and r is the radial distance. The thicker the disk, the 
faster it accretes. 

The outline of the paper is as follows. In § 0, we provide an overview of the novel 
technique for forming evolutionary sequences out of hydrostatic models of thick accretion 
disks. The method relies primarily on global conservation laws and the assumption that the 
viscosity is everywhere positive. In § |^, we describe the hydrostatic models which compose 
the evolutionary sequence. The method for evolving the models in time is presented in § ^. 
Initial conditions are given in § ^ for the evolutionary sequences. In § |^, we present results, 
including discussions of the energy spectrum, the residual mass, and the importance of the 
inner nozzle. Long-term emission properties of thin disks which form after the thick disk 
phase are described in § |^ Conclusions are given in § |. 

2. Method Overview 

Here we investigate what ensues after the formation of a thick disk created from the 
debris of the tidal disruption of a star by a massive black hole. In such an event, it is 
believed that a fraction of the stellar debris will be deposited in an accretion disk around 
the massive black hole. 

The evolution of these thick, radiation-pressure-dominated disks is studied as they 
generically expand, accrete onto the black hole, and transfer angular momentum outwards. 
Often, the thick disk expands, begins to accrete, and after a fraction of the disk accretes 
onto the black hole, the disk becomes thin. The evolution is followed until the disks become 
thin. 

The method is to construct hydrostatic models of thick disks using a small number 
of parameters to specify the structure of the disk. These models are then formed into a 
time sequence in which the global quantities of mass, angular momentum, and energy are 
conserved. Viscosity is required to be positive, but not constant as a function of radius or 
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time. The run of viscosity with radius is set by the time sequence of models, as opposed to 
conventional methods in which the standard viscosity parameter, a, is specified and taken 
as constant throughout the disk. Time varying a has been commonly invoked in the study 
of dwarf-novae (e.g. [Uannizzo, Shatter, fc Wheeler 1988|, [Livio fc Spruit HJUT]) as has radial 



variation in a (ISmak 19841 , [Cannizzo, Wheeler, fc Polidan 1986]) ; this paper's method is 



well suited to this general case. The range of allowable time sequences is investigated to 
determine limits on relevant physical parameters such as the efficiency of conversion of 
gravitational binding energy to radiation and of the amount of disk mass left after the thick 
disk phase. It is also possible to evaluate the importance of initial values to the problem 
such as the pericenter of the disrupted star and the mass of the black hole. Our approach 
is similar to that of Abramowicz, Henderson, and Ghosh (1983) who formed sequences of 
hydrostatic thick disk models, by first assuming a constant viscosity parameter, a, and then 
evolving the disk according to global conservation laws. 

The hydrostatic models for non-accreting disks have four free parameters: the inner 
radius of the disk, a polytropic constant which is related to the ratio of gas pressure 
to radiation pressure, and two additional parameters which specify the run of angular 
momentum with radius. For accreting disks, it is necessary to specify only one parameter 
for the run of angular momentum. The models are described in detail below. 



3. Hydrostatic Models 

For most the models we consider, the timescale for evolution of the disks is much longer 
than the orbital time, so we construct hydrostatic models (see, for example, Paczyhski and 
Wiita 1980, hereafter PW). Neglecting the time dependence in the disk equations is a good 
approximation unless the local viscous term, a, approaches unity. Additionally, there are 
non-hydrostatic effects of accretion at the inner nozzle which we discuss in 



To approximate the effect of the Schwarzschild geometry, we use the gravitational 
potential described in (PW): 

. . _ Ji^, . 2^^. (3) 

R — Rs 

The Keplerian angular momentum in this potential is 

jKep = (GMi?)V2_^. (4) 

n — Its 

The Newtonian potential and angular momentum are —GM/R and (GMRY^'^ . Both the 
Schwarzschild geometry and the PW potential differ from the Newtonian potential in that 
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there is a marginally stable orbit, — SRs, and a marginally bound orbit, i^mb — 2it!s- 
The marginally stable orbit occTirs because angular momentum increases with decreasing 
radius when R < R^^] the binding energy, or the sum of the kinetic and potential energies, 
reaches a minimum at -Rms- The binding energy reaches zero at -Rmb- Keplerian thin disks 
can extend inwards only as far as -Rms, whereas thick disks can reach as far as i?mb because 
of pressure gradients. We make no further attempt to account for other general or special 
relativistic effects, which we beheve do not quahtatively affect any of the paper's results 
and may only slightly affect the quantitative results. (For example, the accretion efficiency 
of a thin disk is 5.7% whereas the PW potential yields an efficiency of 6.3%.) 

The first free parameter of the hydrostatic models is the inner radius, ri, where 
T^i > Rmh- The structure of a thick disk is fully determined by the run of angular momentum 
with radius, j{r), at the surface. In order for a disk to be dynamically stable angular 
momentum must increase with radius, so we require the angular momentum distributions 
to satisfy this constraint. 

For the models studied here, we take a modified power law angular momentum 
distribution as a function with two parameters, Tc and b: 



where the constant rb/rs = 1/(1 + &) ensures that dj/dr > for all r. Examples of angular 
momentum distributions are shown in figure 1. The slope, b, ranges between 0.5, where j{r) 
is Keplerian, and 0.0, where the angular momentum distribution is fiat. The constant. A, is 
chosen so that the angular momentum is equal to Keplerian at a specified central radius, r^. 
For r < Tc, the angular momentum is greater than Keplerian and pressure gradients apply 
a force inwards. For r > Tc, the angular momentum is less than Keplerian and pressure 
gradients apply a force outwards. 

For hydrostatic thick disks which are (slowly) accreting, there is a natural constraint 
on the angular momentum distribution (Abramowicz, Jarosynski, and Sikora 1978). The 
angular momentum at the inner radius, ri, must be Keplerian and ri < 3i?s = -Rmb- For a 
non- accreting thick disk, the inner radius can extend inwards of -Rmb only as long as the 

> jKep{f)- When j(r) falls below jKep{f): orbits become unstable. 

The last free parameter is related to the equation of state of the gas. It is convenient 
to assume a polytropic equation of state. Here we take 



j(r) = yl(rc,6)r* ^'^jKep{r) for r > 
j{r) = j(^b) for r < rb. 



(5) 
(6) 



P = Kp''\ 



(7) 



where K, the polytropic constant, is taken to be the same everywhere in the disk at a given 



instant of time. This approximation is explained in § ^4.3| . The barytropic equation of state 
enforces uniform rotation on cyhnders. 

Equations for hydrostatic equihbrium can be solved as a function of the four parameters, 
Tijc, b, and K. We adopt a convention of lowercase r for cylindrical radial distance and an 
uppercase R for the spherical radial distance. 

In the equatorial plane, the radial force balance equation is 

,2 



1 9P _ 9^ 

p dr dr 



V 

r 



where 



The vertical force balance is 



3 

V = -. 

r 



. ,^,d^ IdP 



0. 



dP 1 /, P 

— = 4:Kp^^^ = 4-, 

P P 



The enthalpy, H(r), can be introduced and evaluated as 

H{r) = 

for our equation of state (Eq. ^. With equation | this can be written, 

r ?^ 

H(r) = '^(r) - '^{vi) + / ^dr. 



(9) 
(10) 



(12) 



The differential equation |10] can be integrated to yield 

ci GM 





■ 1 ■ 




p = 







4 A{R - Rs) 



where Ci is constant on cylinders given by 

1/3, . GM 



ci = 4Kp'o^'{r) 



r-Rs 

and pQ is the density in the equatorial plane. 
The disk surface is given by 



H(r) 



GM 
r-Rs' 



(13) 



(14) 
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With the density and run of angular momentum determined, most physical quantities 
of the disk can be obtained from integrations (either numerical or analytical). The disk 
mass and angular momentum are 



M = / 27rrSdr (16) 
J = r2nrEjdr (17) 



The total disk energy is the sum of potential, kinetic, and internal energies (self-gravity is 
neghgible): 

= r?^dr (18) 

pre fZ^ax GMo 

Epot = / 2 / 27rr- ^dzdr (19) 

= / 27rr(3 - -/3)2 / Pir,z)dzdr, (20) 

Jri 2 Jo 

where j3 is the ratio of gas pressure to total pressure and is directly related to the polytropic 
constant, K. The disk luminosity can be determined if we assume that the disk atmosphere, 
which is radiation pressure dominated, is radiating critically like high-mass stars, so they 
can be described by an Eddington model: 

F,^ = {l-/3)-g,s, (21) 

K 

where k is the opacity per gram and g^s is the effective acceleration (including centrifugal 
force). Integrating over the surface of the disk, the luminosity is (PW) 

Lrad = 2 / da (22) 

Jo Jn 

AtiGMc 



K Jr, 



Gm^zr^ GMzr z{R- Rs)^\ ' ^ ^ 

where a is the element of surface area. The luminosity may be altered by strong winds 
(e.g. U97) as well as by effects in the funnel (e.g. [Nityananda fc Narayan 1982| ), but for the 
purposes of this paper we calculate the luminosity with Eq. |2^. Figure 2 and figure 3 show 
many examples of these hydrostatic models. 



4. Evolution in Time 



Our goal is to form time sequences of hydrostatic models under four physical 
constraints: global mass conservation, global angular momentum conservation, global 
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energy conservation, and positive local viscosity. A computer code was written to perform 
the necessary numerical integrations for mass, angular momentum, energy, etc. Additionally, 
it imposes the physical constraints and forms time sequences of models. The physical 
constraints greatly reduce the degrees of freedom in the four dimensional parameter space 
as described below. 

Mass and angular momentum conservation require that for a non-accreting disk M and 
J are constant between models in a time sequence. For an accreting disk, M and J must 
vary as 

A J = j(ri)AM. (24) 

As discussed above i{r\) = jKep(n), for accreting disks. The mass of the black hole, 
Mbh ^ M, so we take Mbh = const. The code ensures that the constraints are meet to 
approximately one part in 10^ per step which far exceeds the necessary accuracy as most 
evolutionary sequences contain about 100 time steps. The constraints of conserved mass 
and angular momentum are not necessarily met if winds carry much the disk mass away 
(e.g. U97). If this is the case, it will be very difficult to detect tidal disruption events from 
their long-term emission, but the bright flares may be even more optically accessible as the 
winds may reprocess some of the higher energy photons to optical. 

Energy is not conserved in the black hole/disk system because of radiation. However, 
we use energy conservation to determine the amount of energy radiated between two models 
and thereby find the time step. For a non-accreting disk. 

At = ^— , (25) 



where AE is determined from equations |1^ and we impose the physical requirement that 
AE < 0. For an accreting disk, the energy of the system changes not only due to radiation 
but also due to mass loss on to the black hole. In contrast to advection dominated models 
(e.g Narayan & Yi 1995, Chen et al. 1995) the pressure at the inner radius falls to zero, so 
accreting matter has negligible internal energy. Therefore, the energy change between steps 
is 

-GM , j{nf' 



AE = -Erad + AM 



(26) 



The constraint of positive viscosity ensures that angular momentum is transfered 
outwards. Specifically, for angular momentum distributions for which dQ/dr < which 
are of general interest in astrophysics, the torque between two adjacent cylinders is always 
such that angular momentum is fiowing outwards on mass shells. This constraint means 
that the torque between adjacent cylinders, g, is greater than zero, so J — Mj > 0, since 
J = Mj + g. Thus, in contrast to the first three constraints, the viscosity constraint needs 
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to be applied locally. In practice, the constraint is applied at twenty points in the disk 
which are evenly spaced in mass and include the endpoints. 

In geometric terms, the goal is to understand the physically allowable paths that can 
be taken from an initial starting point in the 4-dimensional parameter space of ri, the inner 
radius; Tc, the radius at which the angular momentum is equal to the Keplerian value; b, 
the slope of the angular momentum distribution; and K, the polytropic constant. These 
paths are discussed subsequently. 

4.1. Non- Accreting Evolution 

The physical constraints greatly limit the possible evolutionary sequences. The two 
constraints of mass and angular momentum conservation in the non-accreting case reduce 
the parameter space from four dimensions to a 2-d surface (in 4-d) of models with constant 
mass and angular momentum. The remaining two constraints from energy conservation and 
positive viscosity further limit the possible evolutionary sequences, but in a more complex 
manner. 

For purposes of discussion, let us suppose that the energy loss from radiation is 
negligible so that the disk evolves with constant energy. Then, the allowable parameter 
space contracts from a surface to a curve. Finally, the requirement that viscosity is positive 
gives a direction of evolution along the curve, so there is only one physically allowable 
sequence of models which all lie along the curve. 

Figure 4 shows evolution from a narrow torus of matter to a flattened disk in this limit 
with no radiation loss. Initially, when the angular momentum distribution is fiat, the disk's 
surface has a circular cross section for an infinitely narrow torus as can be analytically 
demonstrated. As the disk evolves, the circular cross section evolves into an eUipse as the 
angular momentum distribution steepens. The relationship between the slope of the angular 
momentum distribution, b, and the ellipticity is derived analytically in Appendix B. The 
numerical code well reproduces the relationship predicted by the analytical approximation. 
Like thin-disks, the non-accreting thick disk spreads inwards and outwards in radius as time 
progresses. 

In the more general and physically applicable case where energy is radiated, the 
evolutionary sequences are not unique; however, they all lie on the same 2-D surface. 
The problem can be simplified because the disk is in general spreading both inwards and 
outwards in radius until it begins to accrete, so the code may can take steps at times when 
Ti is equal to specified values less than the starting radius. Once ri is specified, there is only 
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one degree of freedom left. The range of possible models can be bracketed in a number of 
ways. For instance, at each step one can take the minimum and maximum values of the 
angular momentum slope, and fcmax) which provide solutions that meet the positive 
viscosity and energy constraints. Alternatively, the two limiting cases could be taken to 
be models, denoted AE'min and AE'max in which as little or as much energy is radiated as 
possible between timcsteps. Examples of these two limiting cases are shown in figure 5 and 
figure 6. The two types of bracketing cases are identical for much of the evolution of thick 
disks, which is to say the limiting cases for femin are the same as for AEmm- The case with 
^-E'min is one in which as little energy is radiated as possible; and energy is eventually lost 
through low-efficiency accretion. This hmit is of physical interest, because it is the case 
which determines the minimum accretion efficiency of the tidal disruption events as well 
as determines the minimum amount of mass which remains in the disk after the thick disk 
phase. The AEj^^x limit is one in which the disk radiates as much energy as possible. For 
the models wc arc investigating, the slope, b. increases far enough that the disk eventually 
takes a Keplerian angular momentum distribution and becomes essentially thin. While thin 
disks are not well modeled by the approximations of the code (particularly, the polytropic 
constant, K, is not be the same everywhere in the disk), the late stages of evolution of 
a thin disk around a black hole have been studied by Cannizzo, Lee, & Goodman (1990; 
hereafter CLG), and we can draw on their work to determine the subsequent evolution. 

4.2. Accreting Evolution 

While the disk is accreting, there are only three free parameters as discussed above, 
because the angular momentum must be Keplerian at the inner radius, ri. This constraint 
creates an algebraic relationship between ri, b, and Tc. Eliminating Tc, for instance, reduces 
the parameter space to three dimensions. 

Because the disk is accreting, timesteps are taken at specified values of M less than 
-/^initial in & Similar fashion to how timesteps were taken at given values of ri in the 
non-accreting case. The mass and angular momentum constraints, AMn = M-^-i — Mn 
and AJn = j(ri)AMn, can then be used to reduce the three dimensional space to a 
one-dimensional space. As for the non-accreting case, we are left with the task of bracketing 
the possible evolutionary tracks that are allowed within the constraints of energy loss and 
positive viscosity. 

One possible way to bracket the range of time sequences is to take extreme values of 
Ti at each step. This limit is computationally convenient because r; is always between 2 
and 3 Rs in the accreting case. It is a physically interesting limit because the efficiency and 
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disk structure are strongly tied to the inner radius. In particular, the total efficiency of the 
accretion process determined by the effective potential at r; and is given by 

1 /■*/ if r^i 
e = L{t)dt^ -\ 

where an initial disk of mass, mg, accretes at a rate of dm/dt{t) and radiates at a luminosity, 
L(t) from time ti to time tf. As discussed in PW with reference to steady-state models, 
because the effective potential (enclosed in brackets in Eq. ^) rises from a minimum at 
3i?s to zero at the marginally bound orbit, 2Rs, the more mass lost at small ri, the lower 
the efficiency. The full general relativistic treatment for a non-rotating black hole has the 
same qualitative features and is quantitatively well approximated by the potential of PW 

(e.g. emax,s/emax,PW ~ l-l)- 

4.3. Additional Constraints 

Additional constraints can be placed on the disk evolution to make the evolution more 
physically realistic. We discuss two constraints, one of which is related to the local viscosity, 
and the second to the thermal timescale. 

We may wish to place limits on the local or average value of the viscosity. For instance 
when the strength of the viscosity is parameterized with a, i.e. r] = arj^ax = fctp'^s-^o, 
a should always be less than 1. Additionally, when a approaches unity, neglect of 
the time-dependent terms in the momentum-balance equations may become a poor 
approximation. 

Positive viscosity is enforced from timestep to timestep. The viscosity varies with 
radius and time in contrast to most approximations (e.g. standard a-disk theory). The 
local viscosity can be determined from the relation, J = Mj + g, since J and M (averaged 
over the timestep) can be determined at each radius. The torque can be approximated by 
g ^ —a2'n'r^'Evs2zsurdfl/dr. In this manner, a, can be determined as a function of time and 
radius. At each timestep, additional constraints can be placed on a. For instance, a is 
required to be less than 1, and in some cases, the mass weighted value of a is required to 
lie within a certain range. 

A second constraint can be applied to enforce the the assumption that the disk is 
a constant entropy polytrope. The assumption requires that heat be exchanged rapidly 
within the disk, because at large radii, the energy radiated is larger than the energy locally 
dissipated (and vice- versa for small radii). The energy radiated in a shell (Eqs. [2T[ - P^ goes 



BH 



ri{t)-Rs ' 2n{t) 



+ 



dt 



dt 



(27) 



- 12 - 



as rg^s 



r 



. The energy dissipated per unit volume is approximately 




r— — ahVs ~ T 



,26-4 



(28) 



which is generally a sharply falling function of radius, since b is between and 0.5, and E, 
Vg, and even a, generally decrease with radius for r > r^. Energy could be transported with 
the disk either through circulation patterns which are limited by the sound speed, and which 
can be very efficient as shown by contact binaries (see explanation by Tassoul 1992, and 
references therein) or by radiative diffusion which is limited by the diffusion speed, ~ c/r, 
where r is the optical depth between the points of energy transfer. Radiative diffusion 
may become the dominate heat transfer mechanism at very large radii as the density falls. 
Rather than model the radiative diffusion or a circulation pattern, we consider two cases 
with the numerical code: (1) the two mechanisms together are efficient enough to keep the 
disk in a constant entropy state or (2) the heat transfer is limited at a critical radius where 
the time for heat to be transfered by circulation patterns exceeds the Kelvin-Helmholtz 
time. The transfer time is estimated as 



In the second case, we keep track of the amount of matter outside of the critical radius and 
assume that it collapses to a thin disk on the Kelvin-Helmholtz timescale. Once in a thin 
disk, the matter is ignored for the remainder of the evolution. In practice, there is very 
little difference between the main results (§ ^ in of both cases. 



We run the simulations for three black hole masses, 10^, 10^'^, and lO^M© to cover the 
range of black hole masses about which thick disks are likely to form. Because angular 
momentum is lost when the star disrupts, the mean angular momentum of the initial disk 
is the same as that of the star. If the debris circularizes into a narrow torus, then the torus 
would have radius, 2Rp, where Rp is the pericenter of the star's orbit. We consider disks 
with mean angular momentum corresponding to disruptions at 0.5, 1, and 2 Rt, where 
Rt is the the tidal radius. The debris is taken initially to have constant specific angular 
momentum, as might be expected since the debris does have constant angular momentum 
when the star is disrupted and if shocks efficiently circularize the debris, then energy may 
be more readily transfered than angular momentum. We begin the evolution with a high 
internal energy state, where r[ ^ Tc/2, but the results are not very sensitive to this condition 




(29) 



5. 



Initial Conditions 
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(~ 10% difference in the physical parameters discussed below). Nevertheless, there is some 
uncertainty in the initial angular momentum state of the disk. 

The simulations are most applicable to disruptions around black holes with masses 
below ~ 10"^ Mq. The simulations begin with a debris torus which subsequently accretes at 
Eddington and super-Eddington accretion rates; however, the marginally bound disruption 
debris only gradually returns and circularizes on a timescale t^^^ (eq. ^). For black holes 
above ~ 2 x IO^Mq, the circularization timescale becomes so long that material enters the 
disk at a rate slower than required to fuel the black hole at the Eddington rate (U97). For 
these systems, thick disks cannot form. 

For the lower mass black holes, even if the material is returned at a high enough rate 
to fuel the black hole above the Eddington limit, the accretion rate may be so high that 
much of the material accretes before all of the debris is circularized (U97). In this case, 
which appears to generally occur when a ^ 0.01, much of the debris settles into the disk 
while the disk is accreting. This returning matter cannot be directly accounted for with the 
current simulations. 



we 



We address the problem with two types of runs. First, as described in § 
can impose a restriction on the mass weighted viscosity parameter so that the accretion 
times are not so fast as to drain the disk before the material circularizes. These are the 
runs labeled "a" in Tables 1-3. Second, we can evolve a thick disk with as high viscosity 
as physically possible, even though in reality the disk could not have fully circularized. 
Such sequences may approximate the evolution when material is gradually fed to the disk, 
because such simulations have the same mean angular momentum, and it appears that the 
angular momentum is the most important parameter in the problem. In these runs, the 
viscosity generally varies from near unity at the innermost radius to a 0.1 at large radii. 
These runs are marked "c." 



6. Results 

The results from a variety of evolutionary tracks are shown in Tables 1-3. The evolution 
is followed from the initial conditions given in § ^ to a state in which the disk is no longer 
thick. An appropriate cutoff is when the slope, b becomes larger than 0.48, where 0.5 
corresponds to a Keplerian thin disk. The evolutionary sequences are calculated with a 



variety of combinations of the additional constraints of § [4.3| to investigate their effects on 
the main results. 



As explained in § [4.1|, in one type of physically allowed evolutionary track, the disk 
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cools quickly and forms a thin disk. Thin disk evolution has been studied by CLG. If the 
disk cools quickly (a is small), there would likely be a relatively bright period of 1-3 years 
when the disk shines below the Eddington luminosity, followed by a self-similar phase in 
which the luminosity falls off as ~ The long-term evolution of disks is discussed in § |^ 

The main type of evolutionary track we study is one in which the disk accretes quickly. 
In these tracks, the maximum amount of mass is accreted onto the black hole before the 
disks become thin. The disks can accrete extremely fast, so if the viscosity is large enough 
in the disks for such quick evolution, the accretion may be limited by circularization 
of material which takes at least t^nm (eq. 0). Similarly, the efficiency of the accretion 
can become extremely small in the idealized case where the all of the debris begins in a 
circularized disk and it accretes at the maximum rate ("b" and "c" in tables 1-3). The 
pericenter of the stellar orbit is important in determining how much of the disk can be 
accreted as can be seen by the strong correlation between residual mass and pericenter, as 
is explained in § |6.2| . For the same reason, there is an inverse relationship between the black 
hole mass and the residual mass. 



6.1. Spectra 

The spectra of the disks may be calculated at each time step under the assumption 
that the surface is a thermal. The optical depth is quite large, ~ 10^, even to Thompson 
scattering, so the radiation from the disk surface is most likely thermal. With the the 



local luminosity specified by eqs the local temperature may be determined with the 
Stefan-Boltzmann law. A Planck function is used to describe the local emission, and is 
integrated over the disk surface to determine the final spectrum. If as we believe, the 
scattering cross-section far exceeds the absorption cross-section, the local spectrum could 
be better described with a Wien peak and a modified black body spectrum (e.g. Rybicki 
& Lightman 1979). Given the ad hoc nature of the angular momentum distribution and 
uncertainties in the structure of the inner funnel of the disk (e.g. [Nityananda fc Narayan 



1982|) , such refinements are unwarranted. 



Results of such calculations are presented in U97, along with a description of current 
uncertainties which may affect the energy spectrum. As discussed in U97, the viewing 
angle of the disk is very important because much of the emission is radiated from the inner 
regions of the disk which are not visible to an edge-on viewer. The edge-on disk is ~ 50 
times less luminous than the face-on disk, although because the outer disk is much cooler 
than the inner disk, the difference seen in the optical bands may be significantly less, ~ 5. 
Viewed face-on, the bolometric luminosities are a few times the Eddington limit, but the 
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spectra are very hot, peaking between 50 and 100 eV. The estimated bolometric corrections 
for U and V are ~ 5.5 and 7.5, respectively. 

6.2. Residual Mass 

In most cases a significant fraction of the mass of the thick disk cools and forms a 
thin disk around the black hole, where it accretes on a longer timescale (tables 1-3). The 
main reason why mass cannot be fully accreted in the thick disk is that the specific angular 
momentum of the disk rises as the disk accretes, so the disk must spread to large radii. The 
specific angular momentum of material which accretes onto the black hole is generally much 
smaller than the mean angular momentum of material in the disk. For example, the disk 
formed by a solar-type star with pericentcr equal to the tidal radius which is disrupted by 
a IO^Mq black hole has mean specific angular momentum 10GMbh/c, whereas the angular 
momentum of the accreting material, has specific angular momentum 

(ri/i?s) - 1 c c 

If 50% of the disk accretes, then the average specific angular momentum of the remaining 
material must increase to ~ 16, increasing the mean distance of the disk as well as the 
accretion timescales for the remaining material. More generally, if a fraction, x, of the disk 
accretes with a specific angular momentum, 4GMbh/ c, then 

ji - ./ IG'.Ubii/c 
Ji = --, 

1 — X 

As the specific angular momentum of the disk grows, so does its characteristic radius, r oc jf 
and the accretion timescale, t oc r^/^ oc jf. Therefore, when 90% of the mass is accreted, 
the mass accretion rate has fallen by ~ 1000. Thick disks only exist in steady-state when 
the accretion rate is super-Eddington, so the disks will generally cool to form thin disks 
before accreting a large fraction, e.g. > 90%, of the disk. The exceptions in Table 3 occur 
when the initial angular momentum (and therefore pericenter) of the star is small. In fact, 
when the initial mean angular momentum of the disk approaches that of the Keplerian 
angular momentum at marginally bound orbit, the disk can gain very little specific angular 
momentum as it accretes. The initial mean angular momentum will be this low when the 
encounters are very close with ^ 2.5, which may not be uncommon among disruption 
events around black holes with M ^ IO'^Mq. 



(30) 



(31) 
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6.3. Inner Nozzle 

While hydrostatic solutions have no inner nozzle, physically realistic accreting solutions 
must have an opening at the inner radius. Because of the importance of the inner radius in 
determining the efficiency of the accretion, an understanding of the inner nozzle is vital. In 
general, the structure of the inner nozzle depends on the angular momentum at the inner 
radius, the mass accretion rate, and the polytropic constant, K. 

The problem of an inner nozzle for thick disks has been addressed in a number of ways. 
For instance, Loska (1981), solves vertically integrated equations under the assumptions 
of vertical hydrostatic equilibrium, that zjr^ <^ 1, and that the radial velocity of matter 
is dependent only on the r-coordinate. The critical point for the flow is extremely close 
(< 10-Vi), to the inner radius without the nozzle, for the solutions described by Loska. 
The solutions have no simple scaling relation; however, utilizing the fact that the critical 
point is extremely close to the inner radius justifies a slight simplification of the problem 
which provides scaling laws. 

In appendix B, we apply an approach similar to that of Paczyhski and Sienkiewicz 
(1972) who addressed the Roche lobe nozzle in contact binaries. All of the physical 
quantities calculated in the code may be determined for a disk with a nozzle, with 
knowledge of the height of the nozzle (half-thickness of the disk) which obeys the scaling 
relation (Appendix B): 

oc WK''l\\l^\n - iy^'M^^/\ (32) 
so the dependences on physical parameters, especially mass accretion rates, are weak. 

The disks with nozzles differ only slightly from those without. Even for the fastest 
accreting models, the residual mass (§ |6.2| ) differs only by ~ 10 — 20%. The main difference 
in the models with inner nozzles is that the efficiency of accretion is reduced, in some cases 
by a factor of two because kinetic energy is advected into the black hole. 

7. Long-term Evolution 

After the disk becomes thin, the evolution can no longer be tracked by the thick 
disk evolutionary code. To determine the long-term evolution, we begin with the scaling 
relations derived by CLG for thin disks which, in part, draw upon the self-similar solutions 
of Pringle (1974). The mass accretion rate onto the black hole according to their scaling 
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relations (their eqns 17 and 18) and normalized to their numerical results is 

where Md(0) is the initial mass of the disk, and A is the ratio of mean angular momentum 
of the disk to mean angular momentum of a star with pericenter equal to the tidal radius, 
and is ~ 0.5 — 1.5 initially. If no angular momentum were carried into the black hole, 
y4Md(0) would be constant as the disk accretes. The results of § ^ may be applied to the 
long term evolution by using the final state of the thick accretion disk as the starting 
point for the thin disk similarity solutions. The necessary parameters are the mass and 
mean angular momentum of the residual disk. Using Eq. both the luminosity as a 
function of time and possible transitions to an less-efficient advection dominated regime 
may be evaluated. The evolution of the thin disk luminosity was determined by CLG. 
They characterized the duration of the long-term bolometric luminosity with the time, ti, 
at which Lboi/AfeH = Lq/Mq and found it could be many thousand years. Initially, the 
thin (and thick disks) are extremely hot and radiate only a small fraction, ^ 1%, of their 
luminosity in the optical (see CLG and U97). After 3000 years, CLG find that the peak of 
the broad fiux distribution shifts to 2000 A, so a much larger fraction of light is emitted in 
the optical. Assuming an efficiency of 6% for conversion of mass to energy, as is appropriate 
for thin disks around Schwarzschild black holes, the disk will remain optically bright for 
(CLG) 



where f{t) which accounts for the relative brightness of stars and the disk in the spectral 
band of choice; it is the ratio of the fraction of disk light emitted in the band to the fraction 
of stellar light emitted in the band, and is close to 1 in the optical after 3000 years. 

We may now ask the question of how an early thick disk phase could affect the 
long-term emission by applying Eq. ^ to the thin disks which remain after the thick disk 
phase. The results are shown in Tables 1-3. The timescales for dimming are extremely long. 
If the tidal disruption rate in a typical galactic nuclei is ~ 10~^ yr^^ and the thin disks 
accrete with efficiency 6%, then lower mass black holes would be in a near constant state of 
emission. 

If the accretion rate drops low enough, the accretion disk may become advectively 
dominated. Advection would reduce the efficiency of the accretion process below 6% and 
could move much of the emission to a higher temperature out of the optical band. The 
disk may become advective when the mass accretion rate drops below a critical value. 
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Merit ~ 10 a^LEdd/c^- For a number of advective systems, Merit ~ O.lLEdd/c^ (e.g. Narayan 
1996) The time at which the accretion may become advective according to Eq |3^ is 



. 40 ( iL) (m>] (^) years. (35) 



Below the critical mass accretion rate, the luminosity falls as ~ M^. If the mass accretion 
rate continues to follow the same form as Eq. |3^, then the point at which the bolometric 
luminosity falls below MbrLq/Mq is 

t,« 400 .4 / [-) j j year.. (36) 

Tables 1-3 shown that even in the advection dominated case, the black holes will stay bright 
for an extended period of time. For IO'^Mq black holes, the disk will dim on the timescale 
of ~ 200(0. 1/a) years, whereas for lO^M© black holes, the disk may remain bright 10 times 
as long even when advection dominated. 

The self-similar solution does not account for the loosely bound debris of the initial 
disruption which slowly rains down on the disk as M ~ t^'^ (Rees 1988, Phinney 1989), 
but the return rate of the marginally bound material falls below that of the accretion rate 
from the self-similar solution within the first year. We conclude that the long-term emission 
will be determined by the evolution of the thin disk and not by return of marginally bound 
material in agreement with CLG. 



8. Conclusions 

Only a fraction of the disk mass can be consumed in the initial thick disk phase, which 
is likely to occur for tidal disruptions around 10^ — 10^ Mq black holes. This fraction 
ranges from ~ — 0.99 for disruptions around 10^ to 10'' Mq black holes, and has a typical 
value of 0.5 — 0.9. The amount of the residual mass decreases for larger mass black holes 
and increases with the pericenter of the disrupted star. This residual material will accrete 
in a thin disk over a longer period of time. The existence of the initial thick disk phase 
may reduce the dimming timescale of the disk by a factor of ~ 2 from estimates based on 
thin disks alone. Assuming that an O.SM© initial thick disk, even if the thin disks become 
advection dominated, the black hole mass to light ratio can rise above Mq/Lq = 1 in no 
less than 20 (O.l/a) to 2000 (0.1/a) years following a tidal disruption event, depending on 
the black hole mass and initial orbital of the disrupted star. The long-term emission will 
be most prevalent around lO^M© black holes. If the tidal disruption rates in these galactic 
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nuclei are ~ 10 ^ yr ^, then at least ~ 10% of the nuclei should exhibit the long-term 
UV/optical emission at a level of ~ 10^^ ergs s~^. 

I thank Jeremy Goodman, Bohdan Paczyhski, and David Spergel for many helpful 
suggestions on this work. AU was supported by an NSF graduate fellowship and NSF 
grants AST93-13620 and AST95-30478. 



In this appendix, the evolution of an infinitely narrow torus with constant angular 
momentum is investigated. This example can be evaluated in large part analytically and 
therefore offers a check on the the numerical analysis used elsewhere in the paper. 

A thick disk with constant angular momentum and barytropic equation of state of 
infinitesimal width in a Newtonian potential has equipotentials (including the surface) that 
are circular in cross-section. For example, at the surface of such a circular-cross-section 
disk, 



where the disk is centered on (ri + A) and the torus has radius A = (rg — t\)/2. 

More generally, we show that when the angular momentum is not constant with radius, 
but rather a power law, the equipotentials are ellipses. At the surface of the disk (see 
section ^ 



For simplicity, the calculations are performed in a Newtonian potential, {Rs/ri <^ 1), the 
angular momentum is a power law with radius, (j = Ar''), and the distance between the 
inner and outer radii is 2A. 

Within this framework, the equation for the surface is 



9. Appendix A 
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The amplitude of the angular momentum, A, is such that the angular momentum is 
Keplerian at the center of the disk, r[ + A: 

A = {ri + A)^/'-^ (41) 

The powers of A that are needed for the calculation, expanded to first order in 5 = A/r; are 

A = rl^^-\l + {1/2 -b)6) (42) 
A^ = rl-^\l + {1 - 2b)6) (43) 
^4 _ rf-^\l + {2-Ab)5). (44) 

Writing r as ri + e, where —A < e < A, substituting into Eq. and expanding to second 
order in e/n and A/n yields 

= rf (1 + 2- + 2(1 - 26)^ + 2b^) - (45) 

Replacing e with r — and expanding terms 

= (-l + 26)(rf + r2 + 2rA-2rri-2riA) (46) 
+ {-l + 2b)A^ = {-l + 2b){r -n- Af (47) 



2; 



2 



26 



+ (r - r; - A)2 = A^ (4J 



The shape of the surface (and equipotentials) ranges from a circle when 6 = to a thin disk 
or line when b = 0.5, corresponding to Keplerian rotation. 

As expected, the code utihzed in earlier sections produces ellipsoidal cross sections 
with the correct ellipticity. As a non-accreting disk evolves, the inner-edge moves inwards 
and the outer edge moves outwards. As the disk expands, it flattens. Equation H is a good 
description of disk equipotentials when both Rs/n and (rg — ri)/ri are small. 



10. Appendix B 

In this appendix, we derive a useful scaling relationship for the thickness of the 
nozzle. We operate under the assumptions of Loska (1981), who solved vertically integrated 
equations under the following assumptions: the disk is in vertical hydrostatic equilibrium, 
z/ri <^ 1, and the disk is polytropic. According to Loska, the critical point for the flow is 
extremely close (< lO^^ri), to the inner radius without the nozzle, but the solutions have 
no simple scaling relation. The proximity of the critical point to the inner radius justifies 
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a slight simplification of the problem which provides scaling laws. We follow an approach 
similar to that of Paczyhski and Sienkiewicz (1972) who addressed the Roche lobe nozzle in 
contact binaries. 

In absence of accretion, the disk forms a cusp at n. We expand the effective potential 
around this point, fixing r = r[: 

= 0o[l + a^2] + |^, (49) 
2ri 

where the gravitational potential (PW) is 0o = —GM/{ri — Rs). At the disk surface, the 
effective potential is 

0, = 0o[l + a2;2] + ^, (50) 
where is the half-thickness of the nozzle, and 

a = — ^ r- (51) 

nin-Rs) ^ ' 

In hydrostatic equilibrium, 

dP , 

— = -d(j). (52) 
P 

Integrating eq. ^ at r = r;, with a polytropic equation of state, P = Kp^^^/^, yields 

(n + l)Kp^/^ + = 0s,hyd, (53) 

where (n + l)Kp^/'^ = H{r) is the enthalpy (eq. |lTD, and 0s,hyd is the effective potential 
evaluated at the the hydrostatic disk's surface. When the material has a constant radial 
velocity, but maintains its barytropic equation of state, eq. ^ becomes a Bernoulli equation: 

+ {n + l)Kp^''' + = const (54) 

If the location of the surface were known, the enthalpy in the equatorial plane can be 
determined at r\ (since the value of p in the equatorial plane is known from vertical force 
balance (eqs. |10|,|13|,|14D). Once determined at ri, the enthalpy can be found throughout the 
plane (eq. [1^) from which all other physical quantities can be derived. The value of the 
enthalpy at is 

^^"'^ = " {rl + ziyi^-Rs ^^^^ 

The thickness of the nozzle may be estimated as follows. The flux of matter through 
the nozzle is limited by 

F = Airr-i / pv^dz <^ Airri / pvdz. (56) 
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The maximum value of pv may be determined by differentiating eq. 



vdv + K 



n 



-P 



and using the constraint that at maximum 



n 



d{pv) = = pdv + vdp. 
Combining the last two equations with the Bernoulli equation yields 

The flux is then limited to 

^ max ^" ' 1 



n + 0.5 



K{n + l) 



n+0.5 



n 



\n + 0.5y 

, \ 71+0.5 



K{n + l 
n \ 



dz 



n + 0.5 J \K{n + l)J Jo 

n+0.5 / 



r 2 _ ^2^n+l/2^^ 



n + 0.5 



n 



Kin + l] 



where 



C, = ^3}^-^ + = .429515 for n = 3 



(57) 



(58) 



(59) 



(60) 
(61) 
(62) 

(63) 



2r(2 + n) 

The true matter flux, m is slightly less than -Fmax? but in the case of Roche overflow 
in binaries, the difference is no more than 20%, a negligible amount considering the 
crudeness of the assumptions. Nevertheless, a direct calibration on m/Fmax with numerical 
simulations would be welcomed. Simply equating m and Fmax gives an estimate of the 
nozzle half-thickness 



,n(-0oa)"+^/^ 



l/{2n+2) 



-Lli±i^n + l/2)-V2l 
47r n Ci 



l/(2n+2) 



(64) 



For n = 3, 



1.56084 



n 1/ 



(65) 
(66) 
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Table 1: Mbh = 10^. Thick disk simulations with (1) the pericenter, Rp, in units of 
the tidal radius, Rt, (2) the total time the thick disk was evolved, in seconds, (3) residual 
mass after thick disk phase as a fraction of the initial mass of 0.5 Mq, (4) the efficiency 
of accretion (eq. ^TD, (5) the time, ti, for a thin disk with 6% accretion efficiency to reach 
-^boi/^BH = Lq/Mq assuming self-similar evolution discussed in § |^, (6) the time, tadv, at 
which the thin disk could shift to an advection dominated mode if a ~ 0.1, (7) the time, 
t2, for a thin disk which is advection dominated with a ~ 0.1 to reach Lboi/MeH = Lq/Mq. 
Runs labeled "a" have a restricted viscosity parameter, a (generally <^ 10^), such that the 
disk does not accrete faster than the debris can circularize. Runs labeled "b" allow for the 
outer disk to cool and become thin as described in § H^Sl and allow the maximum viscosity. 



allow the the maximum viscosity (a = 1 and therefore the least efficient 
accretion. In runs labeled "d," the disks cooled and became thin before any accretion took 
place. 
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Table 3: Mbh = 10^. For description, see table 1. 
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Fig. 1. — Angular momentum as a function of radius for a Keplerian disk (thick solid line; 
eq. H) in the pseudo-Newtonian potential of PW, (eq. ^ for non-accreting disks (thin solid 
lines), and for an accreting disk (dashed line). 
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Fig. 2. — A non-accreting thick disk. The top panel shows equipotentials containing 100 % 
(surface), 95 %, and 50 % of the mass in the r/z plane. The bottom panel shows the mass 
distribution as a function of radius. 
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Fig. 3. — An accreting thick disk. The top panel shows equipotentials containing 100 % 
(surface), 95 %, and 50 % of the mass in the r/z plane. The bottom panel shows the mass 
distribution as a function of radius. 
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Fig. 4. — Evolution of a circular tube into a disk using the evolution code. The code well 
reproduces the analytical results of Appendix B. The top panel shows cross-sections of the 
surface at different times. The bottom panel shows the mass distribution as a function of 
radius. 




Fig. 5. — Evolution of a non-accreting disk in the limit that very little energy is radiated. 
The initial model is shown as a thicker line. The different curves correspond to the disk at 
different times as it expands. 
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Fig. 6. — Evolution of a non-accreting disk in the limit that the maximum energy is radiated 
at each step, so that the disk quickly becomes thin. The initial model is shown as a thicker 
line. The different curves correspond to the disk at different times as it expands. 



